TOPMed imputed genomics enhances genomic atlas of the human proteome in brain, cerebrospinal fluid, and plasma.

Comprehensive expression quantitative trait loci studies have been instrumental for understanding tissue-specific gene regulation and pinpointing functional genes for disease-associated loci in a tissue-specific manner. Compared to gene expressions, proteins more directly affect various biological processes, often dysregulated in disease, and are important drug targets. We previously performed and identified tissue-specific protein quantitative trait loci in brain, cerebrospinal fluid, and plasma. We now enhance this work by analyzing more proteins (1,300 versus 1,079) and an almost twofold increase in high quality imputed genetic variants (8.4 million versus 4.4 million) by using TOPMed reference panel. We identified 38 genomic regions associated with 43 proteins in brain, 150 regions associated with 247 proteins in cerebrospinal fluid, and 95 regions associated with 145 proteins in plasma. Compared to our previous study, this study newly identified 12 loci in brain, 30 loci in cerebrospinal fluid, and 22 loci in plasma. Our improved genomic atlas uncovers the genetic control of protein regulation across multiple tissues. These resources are accessible through the Online Neurodegenerative Trait Integrative Multi-Omics Explorer for use by the scientific community.

Comprehensive proteomics of CSF, plasma, and urine identify DDC and other biomarkers of early Parkinson's disease.

Parkinson's disease (PD) starts at the molecular and cellular level long before motor symptoms appear, yet there are no early-stage molecular biomarkers for diagnosis, prognosis prediction, or monitoring therapeutic response. This lack of biomarkers greatly impedes patient care and translational research-L-DOPA remains the standard of care more than 50 years after its introduction. Here, we performed a large-scale, multi-tissue, and multi-platform proteomics study to identify new biomarkers for early diagnosis and disease monitoring in PD. We analyzed 4877 cerebrospinal fluid, blood plasma, and urine samples from participants across seven cohorts using three orthogonal proteomics methods: Olink proximity extension assay, SomaScan aptamer precipitation assay, and liquid chromatography-mass spectrometry proteomics. We discovered that hundreds of proteins were upregulated in the CSF, blood, or urine of PD patients, prodromal PD patients with DAT deficit and REM sleep behavior disorder or anosmia, and non-manifesting genetic carriers of LRRK2 and GBA mutations. We nominate multiple novel hits across our analyses as promising markers of early PD, including DOPA decarboxylase (DDC), also known as L-aromatic acid decarboxylase (AADC), sulfatase-modifying factor 1 (SUMF1), dipeptidyl peptidase 2/7 (DPP7), glutamyl aminopeptidase (ENPEP), WAP four-disulfide core domain 2 (WFDC2), and others. DDC, which catalyzes the final step in dopamine synthesis, particularly stands out as a novel hit with a compelling mechanistic link to PD pathogenesis. DDC is consistently upregulated in the CSF and urine of treatment-naïve PD, prodromal PD, and GBA or LRRK2 carrier participants by all three proteomics methods. We show that CSF DDC levels correlate with clinical symptom severity in treatment-naïve PD patients and can be used to accurately diagnose PD and prodromal PD. This suggests that urine and CSF DDC could be a promising diagnostic and prognostic marker with utility in both clinical care and translational research.

Systematic proteomics in Autosomal dominant Alzheimer's disease reveals decades-early changes of CSF proteins in neuronal death, and immune pathways.

Background: To date, there is no high throughput proteomic study in the context of Autosomal Dominant Alzheimer's disease (ADAD). Here, we aimed to characterize early CSF proteome changes in ADAD and leverage them as potential biomarkers for disease monitoring and therapeutic strategies. Methods: We utilized Somascan® 7K assay to quantify protein levels in the CSF from 291 mutation carriers (MCs) and 185 non-carriers (NCs). We employed a multi-layer regression model to identify proteins with different pseudo-trajectories between MCs and NCs. We replicated the results using publicly available ADAD datasets as well as proteomic data from sporadic Alzheimer's disease (sAD). To biologically contextualize the results, we performed network and pathway enrichment analyses. Machine learning was applied to create and validate predictive models. Findings: We identified 125 proteins with significantly different pseudo-trajectories between MCs and NCs. Twelve proteins showed changes even before the traditional AD biomarkers (Aß42, tau, ptau). These 125 proteins belong to three different modules that are associated with age at onset: 1) early stage module associated with stress response, glutamate metabolism, and mitochondria damage; 2) the middle stage module, enriched in neuronal death and apoptosis; and 3) the presymptomatic stage module was characterized by changes in microglia, and cell-to-cell communication processes, indicating an attempt of rebuilding and establishing new connections to maintain functionality. Machine learning identified a subset of nine proteins that can differentiate MCs from NCs better than traditional AD biomarkers (AUC>0.89). Interpretation: Our findings comprehensively described early proteomic changes associated with ADAD and captured specific biological processes that happen in the early phases of the disease, fifteen to five years before clinical onset. We identified a small subset of proteins with the potentials to become therapy-monitoring biomarkers of ADAD MCs. Funding: Proteomic data generation was supported by NIH: RF1AG044546.

Proteo-genomics of soluble TREM2 in cerebrospinal fluid provides novel insights and identifies novel modulators for Alzheimer's disease.

Triggering receptor expressed on myeloid cells 2 (TREM2) plays a critical role in microglial activation, survival, and apoptosis, as well as in Alzheimer's disease (AD) pathogenesis. We previously reported the MS4A locus as a key modulator for soluble TREM2 (sTREM2) in cerebrospinal fluid (CSF). To identify additional novel genetic modifiers of sTREM2, we performed the largest genome-wide association study (GWAS) and identified four loci for CSF sTREM2 in 3,350 individuals of European ancestry. Through multi-ethnic fine mapping, we identified two independent missense variants (p.M178V in MS4A4A and p.A112T in MS4A6A) that drive the association in MS4A locus and showed an epistatic effect for sTREM2 levels and AD risk. The novel TREM2 locus on chr 6 contains two rare missense variants (rs75932628 p.R47H, P=7.16×10-19; rs142232675 p.D87N, P=2.71×10-10) associated with sTREM2 and AD risk. The third novel locus in the TGFBR2 and RBMS3 gene region (rs73823326, P=3.86×10-9) included a regulatory variant with a microglia-specific chromatin loop for the promoter of TGFBR2. Using cell-based assays we demonstrate that overexpression and knock-down of TGFBR2, but not RBMS3, leads to significant changes of sTREM2. The last novel locus is located on the APOE region (rs11666329, P=2.52×10-8), but we demonstrated that this signal was independent of APOE genotype. This signal colocalized with cis-eQTL of NECTIN2 in the brain cortex and cis-pQTL of NECTIN2 in CSF. Overexpression of NECTIN2 led to an increase of sTREM2 supporting the genetic findings. To our knowledge, this is the largest study to date aimed at identifying genetic modifiers of CSF sTREM2. This study provided novel insights into the MS4A and TREM2 loci, two well-known AD risk genes, and identified TGFBR2 and NECTIN2 as additional modulators involved in TREM2 biology.

Alzheimer's polygenic risk scores are associated with cognitive phenotypes in Down syndrome.

INTRODUCTION: This study aimed to investigate the influence of the overall Alzheimer's disease (AD) genetic architecture on Down syndrome (DS) status, cognitive measures, and cerebrospinal fluid (CSF) biomarkers. METHODS: AD polygenic risk scores (PRS) were tested for association with DS-related traits. RESULTS: The AD risk PRS was associated with disease status in several cohorts of sporadic late- and early-onset and familial late-onset AD, but not in familial early-onset AD or DS. On the other hand, lower DS Mental Status Examination memory scores were associated with higher PRS, independent of intellectual disability and APOE (PRS including APOE, PRSAPOE , p = 2.84 × 10-4 ; PRS excluding APOE, PRSnonAPOE , p = 1.60 × 10-2 ). PRSAPOE exhibited significant associations with Aß42, tTau, pTau, and Aß42/40 ratio in DS. DISCUSSION: These data indicate that the AD genetic architecture influences cognitive and CSF phenotypes in DS adults, supporting common pathways that influence memory decline in both traits. HIGHLIGHTS: Examination of the polygenic risk of AD in DS presented here is the first of its kind. AD PRS influences memory aspects in DS individuals, independently of APOE genotype. These results point to an overlap between the genes and pathways that leads to AD and those that influence dementia and memory decline in the DS population. APOE ε4 is linked to DS cognitive decline, expanding cognitive insights in adults.

Multivariate GWAS of Alzheimer's disease CSF biomarker profiles implies GRIN2D in synaptic functioning.

BACKGROUND: Genome-wide association studies (GWAS) of Alzheimer's disease (AD) have identified several risk loci, but many remain unknown. Cerebrospinal fluid (CSF) biomarkers may aid in gene discovery and we previously demonstrated that six CSF biomarkers (ß-amyloid, total/phosphorylated tau, NfL, YKL-40, and neurogranin) cluster into five principal components (PC), each representing statistically independent biological processes. Here, we aimed to (1) identify common genetic variants associated with these CSF profiles, (2) assess the role of associated variants in AD pathophysiology, and (3) explore potential sex differences. METHODS: We performed GWAS for each of the five biomarker PCs in two multi-center studies (EMIF-AD and ADNI). In total, 973 participants (n = 205 controls, n = 546 mild cognitive impairment, n = 222 AD) were analyzed for 7,433,949 common SNPs and 19,511 protein-coding genes. Structural equation models tested whether biomarker PCs mediate genetic risk effects on AD, and stratified and interaction models probed for sex-specific effects. RESULTS: Five loci showed genome-wide significant association with CSF profiles, two were novel (rs145791381 [inflammation] and GRIN2D [synaptic functioning]) and three were previously described (APOE, TMEM106B, and CHI3L1). Follow-up analyses of the two novel signals in independent datasets only supported the GRIN2D locus, which contains several functionally interesting candidate genes. Mediation tests indicated that variants in APOE are associated with AD status via processes related to amyloid and tau pathology, while markers in TMEM106B and CHI3L1 are associated with AD only via neuronal injury/inflammation. Additionally, seven loci showed sex-specific associations with AD biomarkers. CONCLUSIONS: These results suggest that pathway and sex-specific analyses can improve our understanding of AD genetics and may contribute to precision medicine.

Proteome wide association studies of LRRK2 variants identify novel causal and druggable proteins for Parkinson's disease.

Common and rare variants in the LRRK2 locus are associated with Parkinson's disease (PD) risk, but the downstream effects of these variants on protein levels remain unknown. We performed comprehensive proteogenomic analyses using the largest aptamer-based CSF proteomics study to date (7006 aptamers (6138 unique proteins) in 3107 individuals). The dataset comprised six different and independent cohorts (five using the SomaScan7K (ADNI, DIAN, MAP, Barcelona-1 (Pau), and Fundació ACE (Ruiz)) and the PPMI cohort using the SomaScan5K panel). We identified eleven independent SNPs in the LRRK2 locus associated with the levels of 25 proteins as well as PD risk. Of these, only eleven proteins have been previously associated with PD risk (e.g., GRN or GPNMB). Proteome-wide association study (PWAS) analyses suggested that the levels of ten of those proteins were genetically correlated with PD risk, and seven were validated in the PPMI cohort. Mendelian randomization analyses identified GPNMB, LCT, and CD68 causal for PD and nominate one more (ITGB2). These 25 proteins were enriched for microglia-specific proteins and trafficking pathways (both lysosome and intracellular). This study not only demonstrates that protein phenome-wide association studies (PheWAS) and trans-protein quantitative trail loci (pQTL) analyses are powerful for identifying novel protein interactions in an unbiased manner, but also that LRRK2 is linked with the regulation of PD-associated proteins that are enriched in microglial cells and specific lysosomal pathways.

Unique genetic architecture of CSF and brain metabolites pinpoints the novel targets for the traits of human wellness.

Brain metabolism perturbation can contribute to traits and diseases. We conducted the first large-scale CSF and brain genome-wide association studies, which identified 219 independent associations (59.8% novel) for 144 CSF metabolites and 36 independent associations (55.6% novel) for 34 brain metabolites. Most of the novel signals (97.7% and 70.0% in CSF and brain) were tissue specific. We also integrated MWAS-FUSION approaches with Mendelian Randomization and colocalization to identify causal metabolites for 27 brain and human wellness phenotypes and identified eight metabolites to be causal for eight traits (11 relationships). Low mannose level was causal to bipolar disorder and as dietary supplement it may provide therapeutic benefits. Low galactosylglycerol level was found causal to Parkinson's Disease (PD). Our study expanded the knowledge of MQTL in central nervous system, provided insights into human wellness, and successfully demonstrates the utility of combined statistical approaches to inform interventions.

Proteogenomic analysis of human cerebrospinal fluid identifies neurologically relevant regulation and informs causal proteins for Alzheimer's disease.

The integration of quantitative trait loci (QTL) with disease genome-wide association studies (GWAS) has proven successful at prioritizing candidate genes at disease-associated loci. QTL mapping has mainly been focused on multi-tissue expression QTL or plasma protein QTL (pQTL). Here we generated the largest-to-date cerebrospinal fluid (CSF) pQTL atlas by analyzing 7,028 proteins in 3,107 samples. We identified 3,373 independent study-wide associations for 1,961 proteins, including 2,448 novel pQTLs of which 1,585 are unique to CSF, demonstrating unique genetic regulation of the CSF proteome. In addition to the established chr6p22.2-21.32 HLA region, we identified pleiotropic regions on chr3q28 near OSTN and chr19q13.32 near APOE that were enriched for neuron-specificity and neurological development. We also integrated this pQTL atlas with the latest Alzheimer's disease (AD) GWAS through PWAS, colocalization and Mendelian Randomization and identified 42 putative causal proteins for AD, 15 of which have drugs available. Finally, we developed a proteomics-based risk score for AD that outperforms genetics-based polygenic risk scores. These findings will be instrumental to further understand the biology and identify causal and druggable proteins for brain and neurological traits.

Genetic architecture of plasma Alzheimer disease biomarkers.

Genome-wide association studies (GWAS) of cerebrospinal fluid (CSF) Alzheimer's Disease (AD) biomarker levels have identified novel genes implicated in disease risk, onset and progression. However, lumbar punctures have limited availability and may be perceived as invasive. Blood collection is readily available and well accepted, but it is not clear whether plasma biomarkers will be informative for genetic studies. Here we perform genetic analyses on concentrations of plasma amyloid-ß peptides Aß40 (n = 1,467) and Aß42 (n = 1,484), Aß42/40 (n = 1467) total tau (n = 504), tau phosphorylated (p-tau181; n = 1079) and neurofilament light (NfL; n = 2,058). GWAS and gene-based analysis was used to identify single variant and genes associated with plasma levels. Finally, polygenic risk score and summary statistics were used to investigate overlapping genetic architecture between plasma biomarkers, CSF biomarkers and AD risk. We found a total of six genome-wide significant signals. APOE was associated with plasma Aß42, Aß42/40, tau, p-tau181 and NfL. We proposed 10 candidate functional genes on the basis of 12 single nucleotide polymorphism-biomarker pairs and brain differential gene expression analysis. We found a significant genetic overlap between CSF and plasma biomarkers. We also demonstrate that it is possible to improve the specificity and sensitivity of these biomarkers, when genetic variants regulating protein levels are included in the model. This current study using plasma biomarker levels as quantitative traits can be critical to identification of novel genes that impact AD and more accurate interpretation of plasma biomarker levels.

Large multi-ethnic genetic analyses of amyloid imaging identify new genes for Alzheimer disease.

Amyloid PET imaging has been crucial for detecting the accumulation of amyloid beta (Aß) deposits in the brain and to study Alzheimer's disease (AD). We performed a genome-wide association study on the largest collection of amyloid imaging data (N = 13,409) to date, across multiple ethnicities from multicenter cohorts to identify variants associated with brain amyloidosis and AD risk. We found a strong APOE signal on chr19q.13.32 (top SNP: APOE ɛ4; rs429358; ß = 0.35, SE = 0.01, P = 6.2 × 10-311, MAF = 0.19), driven by APOE ɛ4, and five additional novel associations (APOE ε2/rs7412; rs73052335/rs5117, rs1081105, rs438811, and rs4420638) independent of APOE ɛ4. APOE ɛ4 and ε2 showed race specific effect with stronger association in Non-Hispanic Whites, with the lowest association in Asians. Besides the APOE, we also identified three other genome-wide loci: ABCA7 (rs12151021/chr19p.13.3; ß = 0.07, SE = 0.01, P = 9.2 × 10-09, MAF = 0.32), CR1 (rs6656401/chr1q.32.2; ß = 0.1, SE = 0.02, P = 2.4 × 10-10, MAF = 0.18) and FERMT2 locus (rs117834516/chr14q.22.1; ß = 0.16, SE = 0.03, P = 1.1 × 10-09, MAF = 0.06) that all colocalized with AD risk. Sex-stratified analyses identified two novel female-specific signals on chr5p.14.1 (rs529007143, ß = 0.79, SE = 0.14, P = 1.4 × 10-08, MAF = 0.006, sex-interaction P = 9.8 × 10-07) and chr11p.15.2 (rs192346166, ß = 0.94, SE = 0.17, P = 3.7 × 10-08, MAF = 0.004, sex-interaction P = 1.3 × 10-03). We also demonstrated that the overall genetic architecture of brain amyloidosis overlaps with that of AD, Frontotemporal Dementia, stroke, and brain structure-related complex human traits. Overall, our results have important implications when estimating the individual risk to a population level, as race and sex will needed to be taken into account. This may affect participant selection for future clinical trials and therapies.

Integrative analysis of eQTL and GWAS summary statistics reveals transcriptomic alteration in Alzheimer brains.

BACKGROUND: Large-scale genome-wide association studies have successfully identified many genetic variants significantly associated with Alzheimer's disease (AD), such as rs429358, rs11038106, rs723804, rs13591776, and more. The next key step is to understand the function of these SNPs and the downstream biology through which they exert the effect on the development of AD. However, this remains a challenging task due to the tissue-specific nature of transcriptomic and proteomic data and the limited availability of brain tissue.In this paper, instead of using coupled transcriptomic data, we performed an integrative analysis of existing GWAS findings and expression quantitative trait loci (eQTL) results from AD-related brain regions to estimate the transcriptomic alterations in AD brain. RESULTS: We used summary-based mendelian randomization method along with heterogeneity in dependent instruments method and were able to identify 32 genes with potential altered levels in temporal cortex region. Among these, 10 of them were further validated using real gene expression data collected from temporal cortex region, and 19 SNPs from NECTIN and TOMM40 genes were found associated with multiple temporal cortex imaging phenotype. CONCLUSION: Significant pathways from enriched gene networks included neutrophil degranulation, Cell surface interactions at the vascular wall, and Regulation of TP53 activity which are still relatively under explored in Alzheimer's Disease while also encouraging a necessity to bind further trans-eQTL effects into this integrative analysis.

Gene co-expression changes underlying the functional connectomic alterations in Alzheimer's disease.

BACKGROUND: There is growing evidence indicating that a number of functional connectivity networks are disrupted at each stage of the full clinical Alzheimer's disease spectrum. Such differences are also detectable in cognitive normal (CN) carrying mutations of AD risk genes, suggesting a substantial relationship between genetics and AD-altered functional brain networks. However, direct genetic effect on functional connectivity networks has not been measured. METHODS: Leveraging existing AD functional connectivity studies collected in NeuroSynth, we performed a meta-analysis to identify two sets of brain regions: ones with altered functional connectivity in resting state network and ones without. Then with the brain-wide gene expression data in the Allen Human Brain Atlas, we applied a new biclustering method to identify a set of genes with differential co-expression patterns between these two set of brain regions. RESULTS: Differential co-expression analysis using biclustering method led to a subset of 38 genes which showed distinctive co-expression patterns between AD-related and non AD-related brain regions in default mode network. More specifically, we observed 4 sub-clusters with noticeable co-expression difference, where the difference in correlations is above 0.5 on average. CONCLUSIONS: This work applies a new biclustering method to search for a subset of genes with altered co-expression patterns in AD-related default mode network regions. Compared with traditional differential expression analysis, differential co-expression analysis yielded many more significant hits with extra insights into the wiring mechanism between genes. Particularly, the differential co-expression pattern was observed between two sets of genes, suggesting potential upstream genetic regulators in AD development.